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INTERACTION  OF  PARTICULATE  AND  TRANSFORMATION  TOUGHENING 

by 

John  C.  Amazigo1  and  Bernard  Budiansky 
Division  of  Applied  Sciences 
Harvard  University 
Cambridge,  MA  02138 

ABSTRACT 

A  theoretical  study  is  made  of  the  interaction  between  the  fracture  toughening  effects  of 
crack-bridging  ductile  particles  and  phase-transforming  particles  embedded  in  a  brittle  matrix.  It  is 
found  that  in  certain  parametric  ranges  the  interaction  can  be  synergistic,  with  an  increase  in 
toughness  produced  that  is  greater  than  the  sum  of  the  increases  that  would  be  provided  separately 
by  the  two  types  of  reinforcement.  Quantitative  results  are  provided  for  the  toughening  in  terms  of 
the  individual  toughening  effects  and  coupling  parameters  that  depend  on  the  properties  of  the 
uncoupled  systems. 

/v 

INTRODUCTION 

The  fracture  toughness  of  ceramic  materials  can  be  increased  through  reinforcement  by 
appropriately  chosen  ductile  metal  particles.  A  fairly  general  theoretical  analysis  of  such  particulate 
toughening  in  brittle-matrix  composites,  based  on  the  assumption  that  the  dominant  mechanism  is 
that  of  crack  bridging,  has  been  given  by  Budiansky,  Amazigo,  and  Evans  (1987).  A  different 
micromechanical  toughening  technique  involves  the  incorporation  of  phase-transformable  ceramic 
particles  (notably  Zr02)  into  the  brittle  matrix.  The  associated  transformation  toughening 
phenomenon,  discovered  by  Garvie,  Hannink,  and  Pascoe  (1975),  has  been  studied  theoretically 
(McMeeking  and  Evans,  1982;  Budiansky,  Hutchinson,  and  Lambropoulos,  1983)  on  the 
presumption  that  the  toughening  is  due  to  the  crack-closing  effects  of  dilatant  phase 
transformations  in  the  wakes  of  advancing  cracks.  It  has  been  suggested  (Evans,  1987)  that 
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ductile  bridging  particles  and  transforming  particles  might  interact  synergistically  when  they  are 
both  present,  producing  an  increase  in  toughness  greater  than  the  sum  of  the  effects  they  would 
provide  separately.2  The  outcome  of  an  elementary  analysis  by  Budiansky  (1986)  indicated  that 
the  toughening  ratio  X  =  K/Km  (where  K  is  the  enhanced  toughness  and  Km  is  the  toughness  of 
the  matrix)  could  sometimes  be  approximately  equal  to  the  product  of  the  individual  toughening 
ratios  due  to  ductile  particulates  and  transforming  particles. 

In  this  paper  a  theoretical  study  is  made  of  the  interacting  effects  of  dilatant  transforming 
particles  and  very  ductile  reinforcing  particulates.  The  analysis  will  reveal  conditions  for  the 
validity  of  the  simple  product  rule  for  the  toughening  ratio,  and  will  also  provide  quantitative 
results  for  the  toughening  ratio  when  these  conditions  are  not  met. 


PARTICULATE  TOUGHENING 


The  study  by  Budiansky,  Amazigo,  and  Evans  (1987)  (denoted  by  BAE  henceforth) 
contains  the  following  result  for  the  crack-bridging  effect  (Fig.  1 )  of  ideally-plastic  particles  that 
obey  the  relation  between  particle  stress  cp  and  crack-face  displacement  v  shown  in  Fig.  2.  With 
Kp  defined  as  the  increased  toughness  in  the  presence  of  particles,  the  modified  toughening  ratio 


is  given  by 


Ap  —  Xp/V  (l- c)  — 


KP 


KmVT-c 


Ap  - 


1+ 


2c  ESvf 


-11/2 


^  K„(l-v2) 


(1) 

(2) 


where  c  is  the  volume  concentration  of  particles,  and,  as  shown  in  Fig.  2,  S  is  the  particle  strength, 
and  Vf  is  the  displacement  at  fracture.  The  Young's  modulus  E  and  Poisson’s  ratio  v  of  the 
composite  are  assumed  here,  for  simplicity,  to  be  the  same  as  those  of  the  matrix  material.  This 
last  simplification  was  not  made  in  BAE,  but  is  needed  now  to  avoid  undue  complication  later 
when  transforming  particles  are  introduced  into  the  composite. 


2  Synergistic  effects  have  e /idently  been  observed  (Becher  and  Tiegs,  1987)  in  ceramic  composites  toughened  by  a 
combination  of  whiskers  and  transforming  particles. 
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An  auxiliary  formula  in  BAE  provides  a  relation  between  Ap  and  the  bridge  length  Lp  at 


fracture: 


cS 

Ap  =  1  + 


!  8LP 
7l(l— C) 


Eqs.  (2)  and  (3)  refer  to  the  situation  in  which  the  bridged  crack  shown  in  Fig.  1  propagates  in  a 
steady-state  fashion,  with  simultaneous  occurrence  of  matrix  cracking  and  failure  of  the  particles  at 
the  end  of  the  bridged  zone. 


TRANSFORMATION  TOUGHENING 


In  the  studies  of  transformation  toughening  by  McMeeking  and  Evans  (1982)  and  by 
Budiansky,  Hutchinson,  and  Lambropoulos  (1983)  (hereafter  denoted  by  BHL),  it  was  assumed 
that  during  steady  crack  propagation  martensitic-type  transformations  occur  in  particles  when  they 


are  subjected  to  a  critical  mean  stress  In  the  so-called  supercritical  case  (discussed  in  detail  in 


T 

BHL)  the  full  volumetric  transformation  strain  0p  is  produced  in  each  transforming  particle  in  a 


zone  A  ahead  oi  the  crack  up  as  well  as  in  its  wake,  as  shown  in  Fig.  3.  Along  the  curved  front 
boundary  C  of  this  region,  the  mean  stress  ojj,  is  attained,  and  the  dilatations  thereby  produced 


persist  in  a  wake  of  height  HTon  each  side  of  the  crack.  With  ct  defined  as  the  volume 
concentration  of  transformable  particles,  which  were  assumed  in  BHL  to  have  the  same  elastic 
moduli  as  the  matrix,  the  key  parameter  governing  transformation  toughening  emerges  as 


1+v  Ect0p 


1-V 


Numerical  results  for  the  toughening  ratio  Xj  =  Kf/Km  were  calculated  by  BHL  for  to  in  the  intei  val 


(0,20),  and  were  extended  by  Amazigo  and  Budiansky  (1987)  to  cover  the  full  range  of  to  for 


which  steady-state  crack  growth  is  possible.  Fig.  4  shows  1/Xt  versus  to  up  to  the  critical  "lock¬ 


up"  value  tOc  =  29.99,  first  discovered  by  Rose  (1986). 


The  zone  size  Ht  is  a  function  of  to  in  the  form 


'  f  I  *•  t 


/  Codes 


i  Z  i  Avai  cl .  or 
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„  Km 

0+v)— 


=  g(CO) 


as  shown  in  Fig.  5.  The  results  of  Figs.  4  and  5  combine  to  provide  the  alternative  representation 


,1  sK 

(l+v)— 


=  f(to) 


plotted  in  Fig.  6.  As  shown  in  BHL,  the  toughening  is  related  to  the  zone  size  by 


or,  equivalently. 


Xj  -  [l+2o)g(co)] 


XT=  [l-2cof(Ci))]‘ 


Finally,  we  note  that  the  toughening  ratio  can  be  expressed  in  terms  of  yet  another  zone- 


height  parameter  defined  by 


(1-v) 


*T  2 

- ~  =  (02g(C0) 

K.  1 


Ec,e^ 


For  sufficiently  small  values  of  w,  the  increase  in  toughness  AK  =  Kx-Km  is  given  by 


AK  = 


E^H7c,eTp 


where 


g(0)  =  (47of3)"1  (11) 

COMBINED  TOUGHENING  EFFECTS 

The  results  for  particulate  and  transformation  toughening  are  deficient  in  various  ways.  A 
few  available  measurements  for  particulate  toughening  versus  bridge  length  (BAE)  are  consistent 
with  Eq.(3)  only  if  questionably  high  values  of  particle  strength  S  are  presumed.  Transformation- 
toughening  data  collected  by  Evans  and  Cannon(1986)  show  that  the  increase  in  toughness  appears 
to  be  proportional  to  VHj,  but  the  experimental  values  of  AK  are  substantially  greater  than  those 


predicted  by  Eq.(lO).  Current  thinking  (Evans,1987)  points  to  the  likelihood  that  shear  stress  may 
play  a  role  at  least  as  important  as  mean  stress  in  triggering  the  phase  transformation.  We 
nevertheless  pursue  the  theoretical  analysis  of  the  interaction  of  particulate  and  transformation 
toughening  on  the  basis  of  the  present  models,  with  the  hope  that  final  results  will  be  useful  when 
they  are  interpreted  in  terms  of  experimental  input  values  of  Ap  and  A,j. 

The  analysis  in  Appendix  A  provides  results  for  the  combined  modified  toughening  ratio 

K 


A  = 


Kmvr^c 


(12) 


in  terms  of  Ap,  Xt,  and  a  coupling  parameter 


P  = 


(l+v)cS 


(13) 


■'m 


that  governs  the  interaction  between  particulate  and  transformation  toughening  when  they  occur 
simultaneously  during  steady  crack  growth  (Fig.  7).  Representative  numerical  results  are  shown 
in  Figs.  8-10  for  A  versus  Xj,  for  Ap  =  2,3,4,  respectively.  The  individual  curves  in  each  figure 

are  for  selected  values  of  the  coupling  parameter  in  the  range  [0,°°]- 

The  limiting  results  for  p=oo  and  0  are  of  special  interest.  As  explained  in  Appendix  B, 


we  have 

A  =  ApXy  for  p  — »<*>  (14) 

and 

A  =[  Ap  +  Xj  -1]1/2  for  p  -4  0  (]5> 

Thus  the  previously  anticipated  product  rule  \  =  Xp^yfor  the  combined  toughening  ratio  holds  in 
the  first  limiting  case.  In  this  limit,  as  well  as  for  sufficiently  large  finite  values  of  p,  transforming 
and  bridging  particles  interact  synergistically,  producing  a  larger  increase  in  toughness  than  the 
sum  of  the  increases  that  would  occur  separately.  On  the  other  hand,  for  p  approaching  zero,  the 
combined  increase  can  be  substantially  less  than  cumulative. 

Unfortunately,  in  view  of  the  uncertainties  mentioned  concerning  the  bases  for  the  analysis, 
an  appropriate  choice  for  the  coupling  parameter  p  is  not  easily  made.  However,  the  results  of  the 


analysis  can  be  recast  in  terms  of  an  alternative  coupling  parameter  that  is  much  more  illuminating. 
In  terms  of  the  zone  height  Ht  for  pure  transformation  toughening,  and  the  bridge  length  Lp  for 
pure  particulate  toughening,  the  new  coupling  parameter  is  defined  as 


Ht(1-c) 


The  counterparts  to  Figs.  8-10  are  shown  in  Figs.  11-13,  with  curves  showing  A  for  various 
values  of  T|.  An  appropriate  choice  for  r\  can  now  be  made  on  the  basis  of  observations  of  the 
separate  toughening  phenomena,  bypassing  the  need  to  estimate  dubious  parameters  in  the 
underlying  theory. 

It  is  striking  that  quite  small  values  of  ri  suffice  to  provide  results  close  to  those  for 
t|  =  p  =  oo,  for  which  the  synergistic  product  rule  applies.  This  is  a  happy  result,  because  it  means 
that  synergism  is  not  precluded  despite  the  fact  that  transformation-toughening  zone  heights  tend  to 
be  smaller  than  particulate  bridge  lengths. 

At  fracture,  the  actual  transformation-zone  size  H  and  bridge  length  L  (Fig.  7)  are  generally 
different  from  their  uncoupled  values  Ht  and  Lp.  It  turns  out  (Appendix  C)  that  H  satisfies 

H  A!  -  Aj 

hto-c)  ‘  4-i  (17) 

Hence,  for  given  values  of  Apand  Xj,  the  quantity  H/[Ht(1-c)]  is  an  increasing  function  of  the 


coupling  parameter,  with  (see  Eqs.(l4,15)) 

H  ,  r 

TT-r. — r  =  1  for  T)  =  0 

HT(l-c) 

=  Ap  for  rj  = 

Fig.  14  illustrates  how  H/[Hj(l-c)]  varies  with  "q  for  the  case  Ap=2  and  Ar=3. 

A  simple  general  formula  for  L/Lp  is  not  available,  but  numerical  results  for  L/Lp  vs.  q. 
again  for  Ap=  2,  Xj=  3,  are  shown  in  Fig.  15.  It  is  shown  in  Appendix  C  that  L/Lp  increases 
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L__  + 

Lp  Ap  — 1 


(18) 


at  T)  =  0  to  L/Lp  =  1  for  T]  — »  °o. 

The  essence  of  the  interactive  process  is  the  increase  in  transformation-zone  height  at 
fracture  that  is  produced  by  the  presence  of  bridging  particles.  At  maximum  synergism,  for  which 
the  product  rule  applies,  the  zone  height  is  amplified  by  the  factor  XP. 

CONCLUDING  REMARKS 

The  possibility  of  synergism  between  particulate  and  transformation  toughening  has  been 
demonstrated  on  the  basis  of  simple  models  for  the  individual  toughening  processes  and  their 
interaction.  The  essential  requirement  for  synergistic  interaction  is  that  in  the  uncoupled  situations 
the  transformation  zone  size  not  be  too  small  relative  to  the  particulate  bridging  length.  However, 
this  requirement  is  not  severe;  ratios  of  zone  height  to  bridge  length  of  1/10,  or  even  less,  may 
suffice  to  provide  synergistic  effects. 

It  would  be  desirable  to  corroborate  these  conclusions  by  repeating  the  analysis  on  the  basis 
of  more  realistic  assumptions  concerning  ductile-particle  constitutive  behavior  and  criteria  for  phase 
tranformation. 
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APPENDIX  A 


INTERACTION  ANALYSIS 

We  contemplate  a  steadily  growing  semi-infinite  crack  (Fig.  7)  subjected  far  from  the  crack 
tip  to  conventional  K-field  stresses 


K 


°aP  “  ^=faP(4»  (^22(0)  =  1) 

corresponding  to  the  standard  Muskhelishvili  potential 


(Al) 


(A2) 


The  crack  is  constrained  by  bridging  particles  over  a  length  L,  and,  as  in  BHL,  uniformly 
distributed  plane-strain  dilatation  of  magnitude 


Q=  |(l+v)ct6p 


(A3) 


is  assumed  to  occur  in  the  shaded  region  A.  It  can  be  shown  that  the  elastic  field  associated  with 
the  smeared-out  crack-face  tensions  cS  supplied  by  the  bridging  particles  is  characterized  by  a 
Muskhelishvili  potential  <j>p  satisfying 


pL  i 

Vz-iVL  j 

\  Z  +  2  °g 

^Z+i<L  j 

d<t>p(Z)  cS 
bZ  ~  rt 


and  the  field  due  to  the  transforming  particles  is  described  by 

EQ 


(A4) 


»r-  EQ,  f  log[(Z 
4n(l-v2)Jc 


1/2  +  zi/2)(z1/2  +  zJ/2)]dY0 


(A5) 


outside  the  transformed  region.  The  integral  is  along  the  curved  front  boundary  C  of  the 
transformed  region.  Substitution  of  the  total  potential 
<j>  =  <j>K  +  0P  +  (Jyp  into  the  jump  relation 

V+(X)=^^[4>(X)]-  (A6) 


gives  the  displacement  of  the  upper  crack-face  along  X  <  0,  and  the  mean  stress  outside  A  is  given 


.  _  4<i+v)  r 

_  3  Ke  I  az 


During  steady  crack  growth,  the  displacement  at  X  =  -L  remains  equal  to  its  failure  value 
Vf,  and  the  transformation  criterion  am  =  ocm  must  be  met  as  Z  approaches  C  from  the  exterior  of  A. 


Assertion  of  these  two  conditions  provides  the  relations 


4(l-v2)KVL  4cS(1-v2)L  2(l+v)ct0p  ^Z^+iVT 

vf= - == - = - = - -  Im  log  — — — —  dY0  (As) 

f  EV2^  "E  3tc  I 


<4.2*122*  > 

3  UlrcZ 


l+v^EctOp  i 


2(l+v)cS 


Vi 


+  Im  log 


Vz-iVT 

Vz+iVT 


1— v  J  9k 


1  ,  1  ,Y 
z1/2+ Zq/2  z1/2  +  zJ/2I  0 


The  I-integral  of  BHL,  modified  to  account  for  the  presence  of  the  ductile  particles,  gives 

K2(1-v2)  K^,(1-v2)(1-c)  t 

- 1 - - - g - +  2cSvf+2Homctep  (A10) 

and  another  connection  between  K  and  Km  that  follows  directly  from  a  combination  of  relations  in 


BAE  and  BHL  is 


(All) 


With  the  introduction  of  the  non-dimensional  variables  w  =  v/D,  z  =  x+iy  =  re^  =  Z/D, 


/  =  L/D,  and  h  =  H/D  in  terms  of  the  characteristic  distance 

^  2(l-c)f  (l+v)K„ 


(A12) 


the  governing  equations  (A8-A1 1)  become 
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(A13) 


CO  _  1  1  1  V 

971  Re  Vz  .1/2  1/2  +  1/2  -1/2  dy°  "  1 

^ z  Jr  l  z  +z0  z  +  zo 


(A14) 


.2  ,  4  [  .  coh 

A  =l+_2pwf+T 


(A15) 


(A16) 


where  co,  A,  and  p  are  defined  in  Eqs.  (4),  (12),  and  (13). 

For  prescribed  values  of  A p,  Xt,  and  p,  the  solution  for  A  was  found  as  follows.  The 
value  of  to  is  known  as  a  function  of  Xt  (Fig.  4),  and  wf  may  be  found  from  the  relation 


* 2  _  i  .  8Pwf 

Ap  —  1  -f  — — 
r  371 


(A  3  7) 


given  by  (A16)  for  to  =  0  (which  gives  Eq.(2)  of  the  text).  The  boundary  C  in  the  z-plane  was 
represented  by  r(<p)  in  the  interval  (-0m ,  <j>m)  and  the  expansion 


r 

r(<t> )  sin(4>)  =  £  bn  sin  (n-1/2)^ 


(A18) 


was  assumed.  The  N+3  unknowns  bn  (n=l,2,..N),  4»m  ,  /  ,  and  A  were  found  by  solving 
simultaneously  Eq.  (A  13),  (A  16),  and  the  assertion  of  Eq.  (A  14)  at  4>j=(j<iWN)  ,  j=0,l,2,...N. 


With  h  =  r(4>m)sin(4>ni),  Eq.  (A  15)  provided  a  check  on  the  value  of  A.  Only  a  few  bn’s  were 
needed  for  adequate  accuracy  in  the  preparation  of  Figs.  (8-10). 

The  alternative  coupling  parameter  T)  (Eq.  16)  is  given  in  terms  of  p  by  the  formula 


n 


4p2  ^>T 
(Ap-1)2 


(A19) 


that  follows  from  Eqs.  (3),  (5),  (7)  and  (13). 


APPENDIX  B 

LIMITING  VALUES  OF  TOUGHENING  RATIO 
For  p— »oo  but  Ap  fixed,  we  can  imagine  S-»«>,  L-»0,  with  SVl  and  Svf  bounded.  It  can 
be  verified  that  in  this  limit  the  potential  4>p  defined  by  Eq.  (A4)  vanishes  for  all  Z*0,  but  as  shown 
by  Eq.  (All),  the  crack-face  stresses  nevertheless  continue  to  affect  the  crack-tip  stress-intensity 
factor  Kup=  Km(l-c)1/2.  Perform  the  following  thought  experiment:  in  the  absence  of  bridging 
particles,  apply  K,  and  let  the  phase  transformations  occur.  The  zone  height  H  =  Hj  will  be  given 
by  Eq.  (6).  Now  introduce  the  particles,  and  note  that  since  the  stresses  do  not  change  for  Z*0, 
neither  will  H;  hence  Eqs.  (A  10),  (2),  and  (6)  give 

K2(1-v2)  K^(l-v2)Ap(l-c)  2&>f(co)K2(l-v2) 


It  follows,  via  (8),  that  A2/(Xj)2  =  (Ap)2,  and  so  X=XpXy. 

Next,  for  the  case  p-»0,  we  pretend  that  S-»0,  L— »«>,  again  keeping  SVl  bounded.  In 
this  limit,  the  particle  potential  prescribed  by  Eq.  (A4)  becomes  indistinguishable  from  one  due  to  a 
K-field  (Eq.(A2)).  Withholding  the  ductile  particles,  apply  K.  and  adjust  its  value  to  make  Klip  = 
Km(l-c)1/2'  By  (5),  the  transformation-zone  height  will  assume  the  magnitude 

"  (l+v)Kml2 

H=g(otf  - (1-c) 


Next,  introduce  the  bridging  particles,  and  add  an  equivalent  increment  to  K  to  keep  H  unchanged. 
The  total  K  will  now  satisfy  the  relation 

K2(1-v2)  K2(1-v2)A2p(1-c)  2ci)g(G))K^,(  1  -v2)(  1  -c) 


fill  re 


APPENDIX  C 


TRANSFORMATION-ZONE  HEIGHT  AND  BRIDGE  LENGTH 
From  Eqs.  (5),  (7),  and  the  definition  (A12)  for  D  ,  we  have 


Hj  9n 
D  ~  4(1— c) 


Xj-1 

co 


Using  h  =  H/D  in  Eq.  (A15),  together  with  Eq.  (A17)  gives 

I  =  £  (a2-a2p)  ( 

and  combining  (Cl)  and  (C 2)  provides  the  result  (17). 

To  calculate  L/Lp  from  the  results  of  the  numerical  procedure  of  Appendix  A,  we  write 

L__  L  H  Ht(1-c) 

Lp  H  Hj(l-c)  Lp 


Using  Eqs.  (16-17)  gives  the  formula 


L  T \l(  A2-Ap 


that  was  used  to  prepare  Fig.  15. 

The  limiting  values  of  L/Lp  follow  from  the  observation  that  setting  to  =  0  in  (A  16)  (or 


using  (3))  gives 


AP  =  1  + 


Combining  this  with  Eq.  (A  16)  gives  the  general  result 


2 co  f  dy 

A-l - Re  -= 

/  L  9tc  jcVz 


But  now  recall  the  argument  in  Appendix  B  for  the  limiting  case  p  -»  0,  wherein  field  stresses 
were  unaffected  by  the  introduction  of  particles  and  the  simultaneous  increase  of  K.  This  implies 
that  for  p  — >  0  the  integral  in  (C5)  keeps  the  same  non-dimensional  value  it  would  have  had  in  the 
case  of  pure  transformation  toughening.  Hence,  for  p  — »  0.  and 


Fig.  4 


Transformation  toughening:  reciprocal  of  toughening  ratio  versus  toughening  parameter 
up  to  "lock-up". 


Fig.  5  Dependence  of  transformation  zone-height  on  toughening  parameter 


